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I. INTRODUCTION 

In recent years a number of studies of phase- and 
delay time distributions for waves reflected by a fi- 
nite ID Anderson chain have been published in parallel 
with other studies of phase- J7|-4j2[ and delay time dis- 
tributions |p|,^3|-^5| for continuous disordered chains. In 
particular, an important result which has emerged from 
the latest studies of delay-time distributions |^],[l4|,[l5) for 
different ID systems, is the universality of their asymp- 
totic form, 



P t {t,L -* oo) = T-^e" 
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which was first derived in Refs [9 and [jl3| (see, in par- 
ticular, the footnote 14 in Ref p) and noting that the 
continuum localization length is actually twice the value 
listed in ({J and (see Sect.II.C below)). Here t is the 
time delay, L c the localization length and k the wavenum- 
ber of a reflected wave. Furthermore, the fact that the 
equation (|l|), which coincides with the distribution pre- 
dicted by random matrix theory (RMT) 16 1, has been 
reobtained by Comtet and Texier |l4j for two different 
kinds of disorder in the continuum case and by Ossipov 
et al. ||] for the Anderson model, strongly supports the 
conclusion that the universality of RMT also extends to 
strictly ID random systems where Anderson localization 
dominates. 

On the other hand, the earlier continuous space treat- 
ments [p|-[TT| of phase- and delay time distributions are 
based on an invariant imbedding description Jl7],[l8| of 
reflection and transmission properties of a ID disor- 
dered system. We recall that the invariant imbedding 
method, developed originally by S. Chandrasekhar, leads 
to stochastic differential equations describing the evolu- 
tion of the amplitude reflection (transmission) coefficient 
as a function of length of a disordered sample, viewed as 
a macroscopic continuum Jl7| . The relationship between 
this method and corresponding studies for tight-binding 



chains described by the Anderson model has remained 
somewhat elusive. 

In view of the extensive use which has been made of 
both models of ID chains, it is useful to ex amin e their 
possible connection in more detail. In part II A of Sec- 
tion H we present an exact derivation of the invariant 
imbedding evolution equations from recursion relations 
between amplitude reflection (transmission) coefficients 
of tight-binding Anderson chains of N and N — l atoms 
S, respectively, for weak disorder. 
Next, in part 
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compare available analytical 
and/or numerical results for reflection phase distribu- 
tions for the Anderson- and invariant imbedding 
[p[- |To| , ^2[ models of disordered chains, in order to see to 
what extent they refl ect t he relationship between these 
models established in [EI A . In a similar spirit we briefly 



discuss in part |II Q the effect of discreteness of the Ander- 
son model on the statistical properties of the reflection 
coefficient in a precise way, for weak disorder. 

In the final section [II we discuss the reflection phase- 



and delay-time distributions analytically for a special 
case of a strongly disordered one-dimensional Anderson 
model. Our specific assumption is that the random site 
energies £j follow a doubly peaked distribution composed 
of narrow peaks centered at large positive and negative 
values £i = ±c. This distribution corresponds to a sim- 
ple tractable case of extreme strong disorder for a tight- 
binding chain. The phase- and delay-time distributions 
obtained in this case differ qualitatively from results for 
the more familiar strong disorder case with a wide sym- 
metric distribution of site energies centered at £j = 



II. TRANSPORT IN FINITE CHAINS 
FOR WEAK DISORDER 

A. Derivation of invariant imbedding from Anderson 
model description 

Here we first recall the derivation of Barnes and Luck 
H of a recursion relation for the amplitude reflection 
coefficient of a disordered chain of N atoms described 
by the Anderson model, in terms of the reflection coeffi- 
cient of a chain with one less disordered atom. We then 
extend this discussion for obtaining a similar recursion 
relation for amplitude transmission coefficients of disor- 
dered chains. These Anderson model recursion relations 
are the starting point for deriving the invariant imbed- 
ding stochastic equation for weak disorder. We also ap- 
ply them in the study of reflection phase- and time-delay 
distributions in t he s pecial case of strong disorder con- 
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sidered in Section III. 

The Anderson model for a chain of N dis- 
ordered atomic sites (of spacing a = 1), 
1 < m < N, embedded in an infinite non-disordered 
chain is defined by the tight-binding Schrodinger equa- 
tion 



<Pn+l + <fn-l + e nVn = E<Pr. 



(2) 



where the ip m denote the wavefunction amplitudes at 
sites m and the sites energies e m are random variables 
associated with the disordered sites 1 < m < N. Finally, 
e m = on the two semi-infinite non-disordered chains 
defined by the sites m > N and m < 1 , respectively. 

Let Rn and Tjv denote the reflection and transmission 
amplitudes for an electron wave incident from the right 
with wavenumber —k, and hence energy 

E = 2 cos k , o < k < 7T . (3) 

Thus, by definition, 

<p n = e- ik <- n - N '> + R N e ik( - n ~ N \ forn>N , (4) 

and 

<p n =T N e- ikn , forn< 1 . (5) 
We then define the so-called Riccati ratios 

Y n = , (6) 

which reduce the Schrodinger equation (||) to the two- 
point recursion relation 



Y = E — f — 
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Yn-l 
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where, from (Q), the variable Yn at the last site N of the 
disordered section is related to Rn by 
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By comparing now the equations (^-0) for the chain with 
N disordered sites with the corresponding equations for 
a chain with N — 1 sites obtained by removing the last 
disordered site (ejv = 0), and noting that in both cases 
Yq = e~ lk , it follows that both chains have the same 
Fl, Yi . . . Yjv_i but different Yjy. This immediately leads 
to a recursion relation between Rn and Rn-i, which is 
given by equation (Q) for n — N where for Y/v we insert 
(||), and for Yjv_i the analog of (||) for the chain of iV — 1 
sites since Y)v_i is the same for the chain with N sites 
and for the chain with iV — 1 sites. 
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Rn = 
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where 
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2 sin k 
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We now extend the discussion of Barnes and Luck 
H to obtain an analogous recursion relation for the 
transmission amplitude Tn- This relation is found by 
expressing Tn — e lk tpi, for the chain of length N 
in terms of the Riccati ratios Y n using the identity 

<pi = n^Ti 1 y~(^ + Rn), and Tn-i = e tk <Pi, for the chain 
of length N — 1 in terms of a similar identity. Then, by 
using the fact that the Riccati variables Y\, . . . , Yjv-i are 
the same for the chains of lengths iV and iV— 1 and taking 
the ratio Tjv/Tjv-i we finally get 
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N 



(1 + Rn)Tn-i 
e ik R N -i 



(11) 



By using (|^) , this expression is conveniently rewritten in 
the form 



e 1 N- 
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The equations (^) and (11-a) describe the variation 
of the reflection and transmission coefficients of a disor- 
dered chain when an extra disordered site is added to 
it. The initial conditions for these recursions are obvi- 
oulsy Rq = and Tq = 1. The unitarity relation for the 
reflection and transmission coefficients, 
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may be explicitely checked for a small N e.g. N = 1, 2 , 
For example, for N = 1 one has, from (H) and (11-a), 
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i pi 
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(13) 



which clearly obey (|T^) and illustrate the growth of the 
reflection coefficient and the corresponding decrease of 
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the transmission coefficient of an incident wave, due to 
scattering by a single random site potential. 

In order to derive the invariant imbedding equations 
for reflection and transmission amplitudes from Eqs (^|- 
|i"T| ), we first take the continuum limit. To this end we 
reinstate the lattice constant a and redefine the chain 
length L = Na, the end-site energy e(L) = en and the 
related variable v(L) = v>n, as well as the reflection and 
transmission amplitudes R(L) = Rn and T{L) = Tjy 
and, finally, we define the finite differences R'(L) — 
i [R(L) - R(L - a)} and T'(L) = I [T(L) - T(L - a)}. 
The reinstatement of a leads furthermore to the replace- 
ment of solutions ip n — e lkn for a non-disordered chain 
by (p n = Lp(na) — e lkna and hence to the replacement 
of k by ka in ([|) and in (PHlll). In the continuum 
limit, a — > 0, R'{L) and T'(L) reduce to the correspond- 
ing derivatives dR{L)/dL and dT(L)/dL, so that e.g. 
R(L — a) = R(L) — a dR(L)/dL. Now since the invariant 
imbedding equations are first order differential equations 
it follows that we have to restrict to contributions up 
to order a in the expansion of (^|) and (pT|), where we 
thus approximate e lka by 1 + ika. Using these results the 
continuum limits of (|^) and ( pT| ) become 



a 



dR(L) _ 
dL~ ~ 

2ikaR(L)+w(L)(l+R(L))(l+R(L)+2ikaR(L)) 
l+tu(L)(l+R(L)) 



and 



dT{L) 
dL 



ika + iv{L) 1 + (1 + 2ika)R{L) 



(14) 



dR{L)\ 
(15) 



Next, since the invariant imbedding equations are linear 
in the random potential (proportional to the site poten- 
tial En = e(L) in the present case) we wish to expand 
the r.h.s. of (p^-p^|) to linear order in v^. For this we 
require 



|e(L)| << 2| sin kc 



(16) 



or, within the approximation of (|TJ-^5|), \e{L)\ « 2ka. 
Clearly, this condition is not obeyed near the band edges. 
By expanding the r.h.s. of (p~4jl5h to linear order in e(L) 
these equations may finally be reduced to the form 



ik 



dR(L) 
dL 



= -2k 2 R(L) 
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and 



lk ^H = ^T(L)- £ 44(l + R(L))T(L) , (18) 



dL 
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which coincide with the invariant imbedding stochastic 
equations [ p"7]JT8] ] for the amplitude reflection and trans- 
mission coefficients for a ID chain with a random poten- 
tial 



V(L) = 



e(L) 
2a 2 



(19) 



We emphasize that the validity of Eqs (y_7j|18j) is sub- 
jected to the weak disorder condition ( |l6| ) for wavenum- 
bers k different from the band edge values. On the other 
hand, we recall that previous derivations of these equa- 
tions |l7|jl8| view the disordered chain as an inhomoge- 
neous continuum and assume that the system of length L 
is imbedded invariantly in its extension of length L + AL 
|l7| . Also, these derivations do not involve an explicit 
restriction to a weak random potential. 



B. Analysis of reflection phase distributions in ID 

A variety of detailed numerical and/or analytical re- 
sults for distributions of reflection phases for weakly dis- 
ordered Anderson- or invariant imbed ding m odels have 
been discussed in the literature |p|,[I|,[j| |lC| Ji1| . Here we 
wish to present a more complete comparison of these 
results in the perspec tive of the relation between these 
models established in [II A . Our analysis shows that the 
results for phase-distributions for chains described by the 
Anderson- and invariant imbedding models, respectively, 
are in good qualitative agreement, as expected. 

For the purpose of our discussion we refer to a series of 
figures in the literature where phase distributions for the 
above models have been plotted over a phase interval of 
2ir: we find it useful to classify these figures in tables |(a) 
and |I](b) according to whether they correspond to low- 
6r strong reflection systems, respectively. A weakly re- 
flecting chain of length L (which for the Anderson model 
is identified with the number of atomic sites, N) has a 
localization length L c larger than L while a strongly re- 
flecting one has a localization length smaller than L. Our 
assignment of the numerical results of Refs |p|,^Ul2]| to the 
L << L c or L >> L c categories is based on estimates 
of the localization length using the familiar energy de- 
pendent expression for weak disorder, as given in these 
references. We recall that the curves of fig. 1 of Rcf. p| fol- 
low reflection systems and those of fig. 1 of Ref. pf for 
strong reflection systems are based on explicit analytic 
expressions for the phase distribution while the other fig- 
ures referred to in the tables are obtained from numerical 
calculations 

In the low reflection case the uniform phase distribu- 
tions in figs 2(b) and 2(d), for large kL, of Ref. || for 
the Anderson model are in good agreement with corre- 
sponding nearly uniform distributions of fig. 1 of Ref. || 
for kL = 10 2 and kL — 10 3 , for the invariant imbedding 
model. Similarly, the figs 2(a) and 2(c) of Ref. and 
the fig. 1(b) of Ref. [Q for the Anderson model, for lower 
kL values, agree well with the analytic results of fig. 1 of 
Ref. [pj for the invariant imbedding model, for kL = 5 
and for kL = 10, which indicate a symmetric double 
peak phase distribution in intervals of 2n. On the other 
hand, in the case of strong reflection (for weak disorder) 
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the nearly uniform distributions in the fig. 1 (b) of Rcf . j|] 
and in the fig. 1(a) of Ref. Q for the Anderson model, for 
large kL c are in g ood agreement with the distributions of 
fig. 1 of Ref. jHJ for kL c = 100 and kL c = 200 and with 
the distribution of fig. 6(c) of Ref. Jl| for kL c = 100 
and kL c = 1000, for the invariant imbedding model. 
In the limit kL c << 1 the general expression for the 
stationary phase distribution obtained in has to be 
evaluated numerically. In this limit both the results of 
fig. 1(d) of Ref. [H for the Anderson model and of fig. 6(a) 
of Ref. |l2| for the invariant imbedding model indicate 
that the phase distribution Pe(0) is a symmetric func- 
tion which is strongly peaked about 8 = tt. Furthermore, 
as kL c is increased from zero the peak at 8 = tt shifts pro- 
gressively to larger 8 and a small secondary peak develops 
at 8 < tt, so that the overall shape of the phase distri- 
bution becomes broader, before tending, for kL c » 1, 
to the results discussed above. We note, in particular, 
the similarity between the results of fig. 1 of Ref. (h]] for 
moderately large kL c (kh c — 10 and kL c = 20) and the 
fig. 6(c) of Ref. (H for kL c = 10. 

In conclusion, the above discussion reveals excellent 
overall agreement between reflection phase distributions 
obtained, respectively, form the Anderson model for weak 
disorder and from the invariant imbedding model. This 
is expected to some extent from the equivalence of these 



and 



two models in the continuum limit demonstrated in [I A 



Wc find, in particular, that the discreteness of the Ander- 
son tight-binding lattice does not significantly influence 
the qualitative form of the phase distribution for weak 
disorder. 



C. Reflection coefficient in the Anderson model 



<I*aH 4 > = 



16 sin k 



2N 2 



sin2fciVy 
sin 2k J 



(23) 



In obtaining (|23j) we have expressed averages of prod- 
ucts of four random energies in terms of averages of the 
form (21) using a well-known property of mutually in- 
dependent gaussian variables The corresponding 
results for a rectangular distribution of site energies be- 
tween values — W and W are obtained by replacing Eq by 
W 2 /3. 

The Eq. ( p2[ ) is identified as usual with Ohm's law 
for the Landauer resistance for weak disorder, p ~ 
(\Rn\ 2 ) = 2N/L C , whose expression in terms of the 
localization length follows from (113) and the definition 



\Tn\ 
have 



-N/L e 



L c = 



the limit N « L r . From (\22\i we 



i 2 k 



E = 2 cos k 



(24) 



which coincides with the well-known expression obtained 
by Thoulcss |2(|. Its continuum limit, L c = 8k 2 /eq 
reduces to the localization length L c = 2k 2 /V 2 in a 
continuous gaussian potential with correlation parame- 
ter Vq — —Eq/2 (equatio n |l9|). In fact, in the continuum 
limit (ka — > 0) the Eqs ( p 2| - |23| ) coincide with the results 
obtained previously in an invariant- imbedding treatment 
||. As noted above the continuum limit breaks down 
near the band edges, ka — 0,7r. Finally, we recall that 
the second term on the r.h.s. of ( plj ) describes the effect 
of the inhomogeneity of the reflection phase distribution 



In order to illustrate more precisely the effect of the 
discreteness of the Anderson model for weak disorder in 
a simple typical case, we calculate the first and second 
moments, (|i?Ar| 2 ) and (\Rn\ 4 ) of the reflection coeffi- 
cient. These moments have previously been discussed 
within the invariant imbedding model [^J. In order to 
find moments of |i?/v| 2 for weak disorder we first obtain 
the exact solution for Rn to linear order in the reduced 
site energies (10) by iterating (||) with the initial condi- 
tion i? = 0. This yields 
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Assuming the e m to be independent gaussian variables 
with mean zero and correlation 



F 2 fi 



we obtain successively from (|20|): 
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III. PHASE AND DELAY TIME RANDOMNESS 
FOR STRONG DISORDER 

In this Section we discuss a simple limiting case of 
strongly disordered site energies in the Anderson model. 
We consider a random tight-binding chain where at every 
site the potential takes large values close to a fixed mag- 
nitude, which are positive or negative at random. This 
model corresponds to an extreme case of strong disorder 
since small random site energies are a priori excluded. 
Thus we assume the independently distributed site ener- 
gies £ n to be described by identical doubly peaked dis- 
tributions with peaks centered at E n — ±c, c >> 1, with 
small standard deviation, a << 1. Typical cases of such 
continuous distributions are the doubly-peaked gaussian 
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2crv/27T 
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(25) 



and the doubly peaked rectangular distribution (with 
flat peaks centered at e n — ±c and standard deviation 
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Pe(sn) = [9(e n — c + W) — 8{e n -c-W) 
+ 0(e n + c + W) - 8{e n + c - W)] . 



(26) 



Note that the variances, vare„ = a 2 + c 2 , vare„ = 
W 2 /3 + c 2 ,c » 1, of (H) and @, respectively, char- 
acterize the limit of strong disorder in the present case 
as do large variances in the corresponding single-peak 
distributions (obtained by letting c = 0), vare„ = a 2 , 
with a » 1, and vare„ = W 2 /3, with W/V3 » 1, in 
the usual Anderson model. The advantage of the above 
model of strong disorder is that, due to the fact that typ- 
ical values of the site energies e„ are large, we may use 
perturbation theory in l/e n for studying effects of the 
disorder. 

We now analyse the reflection phase- and delay time 
distributions for the strongly disordered Anderson model 
with the distributions (|25|-p6[) of site energies. We shall 
discuss our analytical expressions in relation to recent 
numerical results for these distributions in the ordinary 
Anderson model for strong disorder B. The probabil- 
ity distribution of the phase of the reflection amplitude, 
Rn — \Rn\e l9n and the corresponding distribution of the 
Wigner delay time, r„ = dO^/dE, will be found from the 
explicit solution of the stochastic recursion relation (||). 
For strong disorder (c >> 1, a « 1 or W/y3 « 1) we 
obtain by successive iteration of (0) 
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Note that, to second order in X/v n , Rn has unit modulus, 
which implies that there is only backscattering to this or- 
der. Furthermore, to order Rn is independent of 
scattering by interior sites n < N and hence independent 
of the length of the disordered chain. In other words, the 
second order expression (|2^) describes reflection in an 
asymptotic stationary regime, L >> L c . 

For simplicity, we restrict to the linear order approxi- 
mation 



Rn — —1 — - — 

iv N 



(28-a) 



for studying the phase- and delay ti me di stribution, 
Pe(ON) and P t (t n ), respectively. From ( 28-a ), ( |Io|) and 
(0) we obtain 



where the distribution of z = 1/gjy isgiven in terms of 
the site energy distributions (|2^) or (g6|), by 
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Finally, from Eqs (|29| - |31| ) we readily obtain the phase- 
and delay time distributions in the form 
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and 
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(33) 



From Eqs (|32|) and (|2j3[|26|) it follows that only small 
phases in the neighbourhood of the peak positions, 9 



± 



2|sinfc| 



have appreciable probability in our model of ex- 
treme strong disorder. We note that the Anderson model 
with an ordinary wide single-peak distribution of site en- 
ergies also leads to a characteristic double peak structure 
for the phase distribution (see the numerical results in 
fig. 1(c) of Ref. Q). However, in contrast to (^), the 
peak positions in the latter case are close to n and their 
separation is of order unity Q . 

Turning to the delay time distribution ( |33"1 ) we observe 
that delay time probabilities are small except near the 
short-time maximum tn ~ ■ — ^ ■ In particular, the ab- 
sence of an asymptotic power law tail, Pt ~ 1/t 2 , as in 
(|l|) , is due to the fact that in the present case an incom- 
ing particle (wave) does not penetrate into the disordered 
sample be yond the outermost site at which it is reflected 
(equation ( 28-a )). 

Ossipov et al. H] have also studied numerically the de- 
lay time distribution for the Anderson model with strong 
disorder i.e. for a wide single-peak distribution of site en- 
ergies. They found a long time tail P T (r) ~ 1/t 2 , t — > 
oo, which is believed to be associated with Azbel reso- 
nances. From this they concluded that the long-time tail 
exists regardless of the strength of disorder. This con- 
clusion is partly disproved, however, by our study of the 
above model of extreme strong disorder. Indeed we find 
that the long-time probability is strongly reduced by an 
exponentially small factor, exp(— c 2 /2a 2 ) « 1, due to 
the fact that the incoming particle is reflected right at 
the edge of the disordered sample. 
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TABLE I. Classification of literature results for Pe(0) pertaining (a) to the low reflection regime, L(N) <C L c , and (b) to 
the strong reflection regime, L(N) ^> L c , in the invariant imbedding model and in the Anderson model for weak disorder. 



(a) L(N) « L c 



Anderson Model 


Reference 


Invariant imbedding 


Reference 


fig. 2(a), N = 10, k = 0, 46tt, L c ~ 6.4 10 b 


2 


fig. 1, kL = 5 


8 


fig. 2(b), N = 50, k = 0, 46tt, L c ~ 2.4 10 7 


2 


fig. 1, kL = 10 


8 


fig. 2(c), N = 10, k = 0, 46tt, L c ~ 2.4 10 3 


2 


fig. l,kL = 100 


8 


fig. 2(d), N = 50, k = 0, 46tt, L c ~ 6.4 10 3 


2 


fig. 1, kL = 1000 


8 


fig. 1(b), k = tt/2, L c ~ 2.4 10 3 


4 




8 


(b) L(N) » L c 


Anderson Model 


Reference 


Invariant imbedding 


Reference 


fig. 1(a), N = 10 3 , k = 271-/3, L c = 290 


2 


fig. l,kL c = 10 


10 


fig. 1(b), N = 10 4 , k = tt/2, L c = 390 


2 


fig. l,kL c = 20 


10 


fig. 1(a), k = ^/ty, L c ~ 60 


4 


fig. l,kL c = 100 


10 


fig. 1(d), k = O.OOlv^F, Lc < 10 -6 


4 


fig. l,kL c = 200 


10 






fig. 6(a), kL c = 0.00001 


12 






fig. 6(a), kL c = 0.0001 


12 






fig. 6(b), kL c = 0.01 


12 






fig. 6(c), kL c = 10 


12 






fig. 6(c), kLc = 100 


12 






fig. 6 (c), kLc = 1000 


12 
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